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A new interpretation of the remarkable X-ray spectrum of 
the symbiotic star CH Cyg 
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ABSTRACT 

We have reanalysed the ASCA X-ray spectrum of the bright symbiotic star CH Cyg, 
which exhibits apparently distinct hard and soft X-ray components. Our analysis 
demonstrates that the soft X-ray emission can be interpreted as scattering of the 
hard X-ray component in a photo-ionised medium surrounding the white dwarf. This 
is in contrast to previous analyses in which the soft X-ray emission was fitted sepa- 
rately and assumed to arise independently of the hard X-ray component. We note the 
striking similarity between the X-ray spectra of CH Cyg and Seyfert 2 galaxies, which 
are also believed to exhibit scattering in a photo-ionised medium. 
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1 INTRODUCTION 

Symbiotic systems contain a red giant star but also ex- 
hibit high-ionisation features in optical/ultraviolet spectra. 
In most cases it is believed that a white dwarf companion is 
accreting from the wind of the giant. 

Symbiotic stars tend also to be X-ray sources and 
iMuerset et alJ lll997l) used R OS AT spectra to separate their 
X-ray emission into three distinct classes: a) supersoft emis- 
sion (below 0.4 keV), jS) emission with a characteristic tem- 
perature of a few 10 6 K (peaking at ~0.8 keV), and 7) harder 
X-ray emission. These classes were interpreted as arising 
from: a) the photosphere of a hot white dwarf, 0) the colli- 
sion of the wind of the giant with a high- velocity wind from 
the hot companion, 7) a possible signature of a neutron star 
companion. One system, CH Cyg, was classed as exhibiting 
bot h a 13 compon e nt an d a "residual 7 component". 

lEzuka et alJ (Il998l) observed CH Cyg with ASCA, 
which was sensitive to harder X-rays than was ROSAT (hard 
limits of lOkeV and 2.5 keV respectively). They found that 
the X-ray spectrum of CH Cyg is made up of apparently 
distinct hard and soft X-ray components. The soft compo- 
nent p eaks at 0.8 keV and the hard component at 5 keV (see 
Fig.Ht. lEzuka et alJ (I1998T) interpreted the hard component 
as heavily-absorbed optically-thin thermal emission of ma- 
terial being accre ted by the white dwarf. In common with 
previous authors llLeahv fc VoUdll995t IMuerset et J ] |l997l) 
they fitted the soft emission as a distinct optically-thin emis- 
sion component. This was interpreted as either coronal emis- 
sion from the giant star or shock heating in the radio/X-ray 
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jets o f CH Cyg. iLeahv fc Volkl (I1995D and IMuerset et all 
il997fl both interpret the soft component as the result of 
shock heating in the collision of stellar winds from the giant 
and white dwarf. 

IWheatlev et al.l (120031) reported a similar two- 
component spectrum in another symbiotic star, 4 Draconis. 
In this case the X-ray spectrum was seen to change be- 
tween the__twa<X5mrMnentform and a single-component 
form. IWheatlev et alJ <l2003l) interpreted both spectra as 
representing emission only from an accreting white dwarf, 
but with a varying amount of partially-ionised absorption. 
Ionisation of the absorbing material removes absorption 
edges entirely from lighter elements (e.g. H, He), and shifts 
absorption edges of heavier elements to higher energies (e.g. 
O), allowing soft photons to leak through. 

In this paper we consider whether similar models can ac- 
count for the apparent two-component nature of the ASCA 
spectrum of CH Cyg. In doing so we draw attention to the 
striking similarity between the X-ray spectra of CH Cyg and 
Seyfert 2 galaxies. 



2 OBSERVATIONS 

CH Cyg was observ ed using the ASCA X-ray observatory 
jTanaka et al.ll994l) on 1994 October 19. ASCA carried four 
X-ray telescopes, of which two were equipped with CCD 
detectors (SISO and SIS1) and two with proportional counter 
detectors (GIS2 and GIS3). The SIS detectors had higher 
spectral resolution and better soft X-ray sensitivity, and in 
this paper we have limited our analysis to data from the SISO 
detector which also had the highest count rate. We have used 
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Table 1. Log of the ASCA SISO observations presented in this 
paper. 

Target seq. no. Date Exposure Count rate 
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Figure 1. ASCA SISO spectrum of CH Cyg. The solid line run- 
ning through the data points represents a simple fit to the spec- 
trum above 2 keV (see text) . The lower panel shows the residuals 
to the fit normalised by the error on each data point. The upper 
line in the top panel shows the unabsorbed model. 



CHCyg 42020000 1994-10-19 19.2 ks 0.516 s" 1 

NGC4507 71029000 1994-02-12 28.3ks 0.126s- 1 

NGC4388 73073000 1995-06-21 31.5 ks 0.032 s" 1 

NGC7582 74026000 1996-11-21 39.7ks 0.011s" 1 



the standard screened data as extracted from the Leicester 
Database and Archive Service (LEDAS) at the University 
of Leicester. This yields an exposure time of 19.2 ks in the 
SISO instrument. 

Data were reduced using FTOOLS software and stan- 
dard techniques as described in the ASCA Data Reduction 
Guide 1 . The SISO source spectrum was extracted from a cir- 
cular region of radius 4 arcmin, and the background was es- 
timated from most of the remainder of the CCD. The mean 
background-subtracted count rate was 0.516 ± 0.006 s~ , 
yielding 9890 ± 120 source counts. A log of observations is 
presented in TableQ 



3 SPECTRAL ANALYSIS 
3.1 Two emission components? 

The SISO spectrum of CH Cyg is pres ented in Fig . DI The 
two X-ray components discussed by lEzuka et alJ 1119981) 
are clearly v isible either side of a minimum at 2 keV. 
lEzuka et alJ 1^998) fitted the spectrum with a complex 
multi-component model consisting of three separate plasma 
emission components (two for the soft emission and one for 
the hard emission) and three separate absorption compo- 
nents (one for the soft emission and two for the hard emis- 
sion). 

The solid line in Fig.0 shows the results of a fit to the 
spectrum only above 2ke V. We used the mekal optically- 
thin t hermal plasma mode l jMewe et alll985lll98 6; Kaastra 
Il992t iLiedahl et all Il995l) absorbed by a neutral cosmic- 
abundance absorber. The best-fitting mekal temperature 
was 7.7 keV (characteristic of accreting white dwarfs) and 
th e absorber co l umn d ensity was 1.6 x 10 23 cm -2 . As noted 
bv lEzuka et alJ il998t) . this simple absorption model does 
not result in a good fit to the hard component (we found 
a reduc ed y 2 of 1.9 with 28 2 degrees of freedom). This mo- 
tivated lEzuka et al.l lll998l) to add a second, complex ab- 
sorption component to their hard X-ray model (they used 
a partial-covering model). They did not, however, consider 
the possibility that the soft component itself could arise as 
a consequence of photons leaking through or around this 
complex absorber. The upper line in Fig.0 shows the unab- 
sorbed model spectrum. It can be seen that the implied soft 
X-ray flux from the hard component exceeds the detected 
flux by an order of magnitude. Even a small soft X-ray leak 
would easily account for the measured soft X-ray emission, 
avoiding the need for a second emission component. 
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3.2 A photo-ionised absorber? 

IWheatlev et alJ <l2003l) successfully fitted the two- 
component X-ray spectrum of the symbiotic star 4 Draconis 
with a single emission component modified by a photo- 
ionised absorber. Ionisation of the absorbing medium 
removes absorption edges entirely from lighter elements 
(e.g. H, He) and shifts absorption edges of heavier elements 
to higher energies (e.g. O). The result is that soft X-rays 
leak through the absorber, giving rise to what can appear 
to be a separate emission component. In this section we 
investigate whether this same model might explain the 
two-component X-ray spectrum of CH Cyg. 

In order to consider the possible effects of ionisation of 
the absorbing me dium in CH Cyg we employed the absor i 
model in xspec iDone et alJ Il992t IZdziarski et alJ Il995l) . 
This is a simple single-zone photo-ionisation model which as- 
sumes a power-law input spectrum. For this paper we modi- 
fied the absori source code to use a bremsstrahlung spectrum 
as the ionising source, and we fixed the bremsstrahlung tem- 
perature to the mekal temperature of the hard component. 
The other parameters of the absori model are the column 
density, Nn, the temperature of the absorber, T a b s , and the 
ionisation parameter £ = L/nR 2 , where L is the integrated 
source luminosity between 5eV and 300 keV, n is the den- 
sity of the material, and R is the distance of the material 
from the illuminating source. In addition to the absori model 
we included a neutral absorption component to account for 
interstellar absorption. 

Fitting with the absori model we found we needed high 
values of £ (~ 10 2 — 10 3 ) in order to allow the soft X-ray leak 
to extend as hard as 2 keV. In contrast, the X-ray minimum 
between the t wo components occurre d at around 1 keV in 
4 Draconis, and Wheatle y et al.l (|2003) found acceptable fits 
with £ ~ 6. 

The best fit with the absori model is plotted in Fig.|2] 
together with the model spectrum. It can be seen that the 
ionised absorber results in a better fit to the hard X-ray 
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Figure 2. ASCA SISO spectrum of CH Cyg fitted with the absori 
ionised absorption model (top panel). The middle panel shows the 
residuals to the fit normalised by the error on each data point. 
The bottom panel shows the model spectrum. 



emission between 2-5 keV than did the neutral absorber 
(Fig.0. It also allows a sufficiently strong soft- X-ray leak 
to account reasonably well for the soft X-ray emission (al- 
though the fit is not statistically acceptable with a re- 
duced x 2 of 2.3 with 348 de grees of freedom). The best- 
fitting column density for the ionised absorber is Ah = 
6.2 x 10 23 cm -2 . The ionisation parameter was found to be 
anti-correlated with the temperature of the absorber, and 
equally good fits were found for £ in the range 55-1000 (for 
Tabs = 10 6 — 10 4 K) . The best-fitting mekal temperature was 
7.1 keV and the neutral absorber had a column density of 
JV H = 2.2 x 10 21 cm -2 . 

Inspection of Fig.[5] shows that the strongest residuals 
occur at 0.9-1.0 keV, coincident with the strong OVIII K- 
shell edge in the absori model. Residuals are also apparent 
around 6.4 keV, the energy of the K a line of neutral iron. 

Adding a narrow emission line fixed at 6.4 keV entirely 
removed the residuals at this energy and resulted in an im- 
provement in reduced \ 2 to 1.94 (with 347 degrees of free- 
dom). The line has an equivalent width of 220 eV, which 
is high for X-ray r eflection from a cold sl ab such as the 
white dwarf surface JCeorge fc Fabianlll99 ll) . This indicates 
either partial ionisation of the reflecting medium or fluores- 
cence from an additional location, perhaps an accretion disc 
and/or the wind of the red giant. 

Removing the residuals in the 1 keV region is more dif- 
ficult. Adding a narrow line also in the 1 keV region resulted 
in a considerable improvement in reduced \ 2 to 1.22 (with 



345 degrees of freedom) with a best-fitting line energy of 
0.95 keV. This line might represent recombination emission 
from oxygen, and indeed the fit statistic was further im- 
proved to 1.12 (345 d.o.f.) by replacing the narrow line with 
a radiative recombination continuum (RRC) with edge en- 
ergy fixed at the OVIII K-shell energy (0.8714 keV). How- 
ever, this fit relies on a high electron temperature (0.1 keV) 
to make the RRC sufficiently broad to cancel out the OVIII 
absorption edge. The perfect cancelling of these two features 
leads us to believe that neither are really present in the spec- 
trum of CH Cyg. 



3.2.1 Fluxes 

The absorbed 0.6-10 keV flux of the best fit to the ionised 
absorber model was 6.3 x 10 -11 erg cm -2 s _1 and the im- 
plied unabsorbed flux of the mekal model was 2.0 x 
10~ 10 erg cm -2 s _1 corresponding to a bolome tric luminos- 
ity of 3 x 10 33 erg s" 1 at a dis tance of 270±65 pc JViotti et alJ 
Il997t [Ferryman et al.lll997l) . This corresponds to an accre- 
tion rate onto a white dwarf of around 3 x 10 16 gs _1 or 
5 x 10 -10 Mq yr _1 (by considering the gravitational poten- 
tial energy of material falling from infinity). This is a fairly 
robust estimate since it depends mainly on the temperature 
and normalisation of the mekal component. It does not rely 
on the ionised absorber interpretation. 



3.2.2 Implied geometry 

As stated above, the ionisation parameter is defined as £ = 
L/nR 2 . The column density is defined as Ah = nl where 
I is the path length through the absorber. The ionisation 
parameter can be written, therefore, also as £ = Ll/NuR 2 . 
The path length through the absorber (I) cannot be much 
larger than the distance of the absorber from the ionising 
source (R), and so there is a maximum allowed distance of 
the absorber from the ionising source is given by i? max ~ 
L/Nu^. Taking the best fitting values for L and Ah and the 
allowed range of £, we find ii max ~ 5 x 10 6 — 9 x 10 7 cm which 
is much smaller than the size of a white dwarf (typically 
~ 7 x 10 s cm). Even allowing for extra X-ray luminosity 
that might be hidden by the absorber, an ionised absorber 
in CH Cyg would have to be located very close to the white 
dwarf. 



3.3 A photo-ionised scattering medium? 

As well as soft X-rays leaking through the absorber, it is also 
possible for soft X-rays to leak around the absorber. This 
is possible if the absorber is not isotropically distributed 
around the white dwarf and instead absorbs preferentially 
along our line of sight. Soft photons emitted in our direction 
would then be absorbed, but soft photons emitted in other 
directions would escape further from the white dwarf where 
they could be scattered into our line of sight in a lower- 
density photo-ionised medium. 

This arrangement is essentially the same as is be- 
lieved to exist in the nuclei of Seyfert 2 galaxies, where the 
anisotropic absorption is attributed to a highly inclined ac- 
cretion disc or molecular torus. The ASCA spectra of Seyfert 
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Figure 3. ASCA SISO spectrum of CH Cyg fitted with xstar 
photo-ionised scattering model (top panel). The middle panel 
shows the residuals to the fit normalised by the error on each 
data point. The bottom panel shows the model spectrum. 

2 galaxies are actually remarkably similar to the spectrum 
of CH Cyg, and we discuss this further in Sect .[I] 

In order to consider the possible effects of photo-ionised 
scattering in CH Cyg we have calcula ted model grids us- 
ing the XSTAR photo-ionisation code l|K allman fc McCravl 
119821: iKallman fc Bautistal EoOll : iTurner et al.lll99rl . This 
code calculates the absorption and emission spectra of a 
spherical shell of material surrounding a source of ionising 
photons. We used the emission spectrum of a model grid 
covering an ionisation parameter (£) range of 10~ 4 — 10 4 
(same definition as the absori model above) and a column 
density range of 10 19 - 10 23 cm . The ionising spectrum 
was chosen to be an 8 keV bremsstrahlung in order to match 
approximately the emission of the white dwarf in CH Cyg. 
This model grid (grid23scatt) is available from the xstar 
webpage 2 . 

We fixed the normalisation of the scattered spectrum to 
that appropriate to the unabsorbed flux of the fitted mekal 
emission component (a value of 4.1 x 10 -4 ). We replaced 
the ionised absorber model with a neutral absorber, but in- 
cluded a second partial-covering absorber in order to account 
for the complexity apparent in Fig.Q We also included a 
narrow emission line fixed at 6.4 keV and a third neutral 
absorber representing interstellar absorption. 

Fitting with this model we found that the photo-ionised 
scattering region readily reproduced the strength and overall 

2 http:/ /heasarc. gsfc.nasa.gov/docs/software/xstar /xstar. html 



flux distribution of the soft component in CH Cyg. We also 
found a statistically acceptable fit, with a reduced \ 2 of 1.09 
(338 d.o.f.). This best fit is plotted in Fig.El 

The strongest residuals in our best-fit scattering model 
occur between 1-2 keV. This corresponds to the energies of 
strong line emission in the mekal model, mainly due to iron 
L-shell, that were not included in our bremsstrahlung ionis- 
ing spectrum. This missing input flux may account for these 
residuals. Alternatively, we have found that relaxing the con- 
straint on the value of the normalisation of the scattering 
component allows a better fit to this region. In effect this 
allows for extra ionising flux that is hidden by the main ab- 
sorber, which is plausible because accretion by white dwarfs 
is known t o result in emission at a wide range of temper- 
atures (e.g. iBaskill et al1l2005h . Finally, these residuals can 
also be reduced by allowing the relative elemental abun- 
dances to vary in the XSTAR model. However, since individ- 
ual emission features are not resolved, none of these possible 
explanations are unique, and we did not pursue them in de- 
tail. We were satisfied that the scattering model is capable 
of producing the strength and overall shape of the soft com- 
ponent in CH Cyg. 



3. 3. 1 Fitted parameters 

The fitted parameters of the scattering region were Nh = 
(5.5±o;i) x 10 20 cm" 2 and £ = 760±Iao- The parameters for 
the mekal emission component were essentially unchanged, 
with kT = 8.3 ± 0.4 keV and an unabsorbed 0.6-10 keV flux 
of (2.29T0.15) x 10" 10 erg s" 1 cm" 2 . The column densities of 
the main absorber are (9.7T0.7) x 10 22 cm -2 and (2.9±0.5) x 
10 23 cm -2 . The former was applied to all the flux from the 
hard component, the latter had a partial covering fraction of 
66 ±4 per cent. The column density of the absorber applied 
to the whole spectrum (including the scattered emission) 
was (2.6 ± 0.5) x 10 21 cm" 2 . The equivalent width of the 
6.4 keV line was 200 ±20 eV. All errors represent 68 per cent 
confidence intervals and account only for statistical errors. 



3.3.2 Implied geometry 

In Sect . I3.2!2*l we showed that our fits to the ionised absorber 
model placed tight constraints on the geometry of the ab- 
sorber. Our scattering model relaxes these constraints be- 
cause the absorption of the hard component is decoupled 
from the source of the soft component. Since the ionisa- 
tion parameter of the absorber is consistent with zero, the 
absorber can be placed at any distance from the ionising 
source. However, the scattering model does require the ab- 
sorption to be anisotropic around the white dwarf, and so 
the absorber is still likely to be closely associated with that 
object. 

The scattering region, in contrast, does have a well de- 
fined ionisation parameter in our model, and so we can 
apply the same argument used in Sect . 13. 2^21 to constrain 
its distance from the ionising source. In this case J? max ~ 
7 x 10 9 cm, which is around an order of magnitude larger 
than the white dwarf. 
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4 DISCUSSION AND CONCLUSIONS 

The spectral analysis of Sect.|^] shows that the soft X-ray 
component in CH Cyg can be interpreted as a reprocessing of 
the hard X-ray component. The X-ray spectrum is explained 
without the need for a second source of X-rays in the system, 
such as from colliding winds. 

In Sect . 13.21 we investiga t ed th e ionised absorber 
model used by Whcatle v et alJ lj2003f) to explain the two- 
component X-ray spectrum of 4 Draconis. This model re- 
produced the strength of the observed soft component, but 
it predicted a strong OVIII absorption edge that is not de- 
tected. The required combination of high column density 
and high ionisation parameter also implied that the absorber 
must be much smaller than the white dwarf. Overall we do 
not regard this model as an adequate explanation of the soft 
X-ray component of CH Cyg. 

In Sect . 13.31 we considered an alternative explanation for 
the soft X-ray component in CH Cyg: scattering of the hard 
component in a photo-ionised medium. In this interpreta- 
tion the absorber is not isotropically distributed around the 
white dwarf and instead absorbs preferentially along our line 
of sight. Soft photons emitted in our direction are absorbed, 
but soft photons emitted in other directions escape further 
from the white dwarf where they can be scattered into our 
line of sight in a lower-density photo-ionised medium. This 
scattering model provided a much more satisfactory expla- 
nation for the soft component in CH Cyg. It readily re- 
produced the strength and shape of the soft component. It 
also implied a distance to the scattering region that is much 
larger than the white dwarf. 

As noted in Sect . 13.31 the ASCA spectrum of CH 
Cyg bears a striking resemblance to those of Seyfert 2 
active galaxies. Examples inclu d e NGC 4388, NGC 7582 
and NGC 4507 ilwasawa et alJ Il997t IXue et al.1 1 19981 
IComastri et alJl998h . In Fig.Hlwe plot the ASCA SISO spec- 
tra of these systems, together with the spectrum of CH Cyg. 
The Seyfert 2 spectra have been extracted using the same 
method as described for CH Cyg in Sect.|5J and details of 
the observations are given in TableQ 

In Seyfert 2 galaxies our line of sight to the central black 
hole is believed to pass through a highly inclined accretion 
disc, or molecular torus, that provides the strong X-ray ab- 
sorption apparent in the ASCA spectra (in some systems 
this region is Compton thick and blocks even hard X-rays). 
Soft X-ray components very much like that of CH Cyg are 
also observed (Fig.[IJl, and this has generally been inter- 
preted as scattering in photo-ionised cones above and below 
the accretion discs/tori. This interpretation was strength- 
ened by the apparent detection of the same ionised material 
in absorption in Seyfert 1 galaxies. The geometry described 
here is illust rated by the schematic diagrams in F ig. 7 of 
lAwakil dl997l) and Fig. 6 of lKinkhabwala et alJ (l2002f) . and is 
essentially the same interpretation as we favour for CH Cyg. 
There has been some debate over whether the soft X-ray 
components of Seyfert 2 galaxies can instead be dominated 
by collisionally ionised gas, perhaps from shock heating in 
outflows from the central region or from nuclear star-forming 
regions (similar interpretations to those proposed by previ- 
ous authors for CH Cyg). However, the scattering model 
has now been confirmed in the Seyfert galaxies Mrk3 and 
NGC 1068 using high-resolution X-ray spectra from Chan- 
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Figure 4. ASCA SISO spectrum of CH Cyg compared with those 
of three Seyfert 2 galaxies. 



dra and XMM-Newton JSako et al.l2000tlKinkhabwala et alJ 
12001 iBrinkman et al]l2002l) . 

The success of the scattering model in explaining the 
X-ray spectra of Seyfert 2 galaxies gives us increased con- 
fidence that this model is a likely explanation for the sim- 
ilar X-ray spectrum of CH Cyg. In Seyfert 2 galaxies the 
source of absorption is believed to be an edge-on accretion 
disc or torus. A highly inclined accretion disc may also ac- 
count for the observed absorption in CH Cyg. The presence 
of a disc is supported by the dete ction of jets llTavlor et alJ 
119861: iGallowav fc Sokoloskill2004l) and there is further evi- 
dence that CH Cy g may be viewed at high inclination (e.g. 
ISkopal et allll996h . 

The scattering interpretation of the X-ray spectrum of 
CH Cyg is attractive because it requires only one source of 
X-ray s: an accreting white dwarf. In contrast, lEzuka et alJ 
( 1998) invoke a second X-ray source, and still require a com- 
plex treatment of absorption and emission (three compo- 
nents of each). Several mechanisms have been invoked in 
order to explain this second source of X-rays, including col- 
liding winds, jet shocks and coronal emission from the giant 
star llLeahv fc VolkllT995l : iMuerset et"ai]|l997l: lEzuka et all 
1998). A consequence of our analysis is that these sources are 
no longer required to explain the X-ray spectrum of CH Cyg. 
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A definitive test of this model would be provided by high- 
resolution X-ray grating spectra of CH Cyg. 
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